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Global topological defects described by real scalar field in (3,1) dimensions coupled to gravity are 
analyzed. We consider a class of scalar potentials with explicit dependence with distance, evading 
Derrick’s theorem and leading to defects with spherical symmetry. The analysis shows that the 
defects have finite energy on flat space, contrary to the observed for the global monopole. With 
the aim to study the gravitational field produced by such defects, after an Ansatz for the static 
metric with spherical symmetry, we obtain the coupled system of Einstein and field equations. On 
the Newtonian approximation, we numerically find that the defects have a repulsive gravitational 
field. This field is like one generated by a negative mass distributed on a spherical shell. In the 
weak gravity regime a relation between the Newtonian potential and one of the metric coefficients 
is obtained. The numerical analysis in this regime leads to a spacetime with a deficit solid angle on 
the core of the defect. 


I.INTRODUCTION 

Topological defects are solutions of the classical field 
equations, which are similar to particles. In quantum 
theory, they correspond to extended particles, composed 
of the elementary particles in each particular modelQ. 
Topological defects are found very frequently in con¬ 
densed matter physics. For instance, the vortexQ, the 
simplest soliton in gauge theory with scalars occur in 
type-2 superconductors. In the cosmological context we 
have cosmic strings, and some models use them as seeds 
for the formation of structures in the early universe^. In 
this way topological defects are of great interest in high 
energy physics |3j. 

Topological defects can be global or local. The global 
defects arise in models with a continuum but global sym¬ 
metry, and have been investigated for instance in0, |b|. 
The local defects usually arise in models which develop 
spontaneous symmetry breaking of some local symmetry. 

As is well known, however, there is a theorem which 
puts a limit on the dimensionality of systems constructed 
only with scalar fields, in order to ensure the presence 
of topological defects^]. Along the years, several dis¬ 
tinct routes have been constructed to evade this theo¬ 
rem, such as: i) to consider constraints on the scalar 
fields, as in the nonlinear 0(3) model|7j; ii) to add gauge 
fields, as with the ’tHooft-Polyakov monopole 0 |9|; iii) 
to include higher-order derivatives on the Lagrangian|ld: 
iv) to introduce time-dependent and non-dissipative so¬ 
lutions (see ref. U3 pp.141-150 and refflil) ; v) to suppose 
non-local solutions, such as global monopoles - this means 
models with a divergent energy that can be limited by a 
cutoff in a cosmological context due to the presence of 
another defect Jj, ,ij. 

Recently, another possibility was considered on search¬ 
ing for defects described by only one scalar field in (3,1) 
dimensions. For this were introduced models described 


by a dimensionless Lagrangian density with an explicit 
dependence with dista,nce|_K|: 

^ = (i) 

where f(x 2 ) is in principle an arbitrary function of the 
dimensionless coordinate <f> a real dimensionless scalar 
field and W = W(<j)) a smooth function of 

The equation of motion for static field is, for a specific 
choice of /(a; 2 ), 


di&$ = v 2 ^ 


1 8Wd 2 W 
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with f = (x^+x^+x 2 ) 1 ^ 2 a radial dimensionless quantity. 
Following R,ef.|13| we look for solutions depending only 
on the radial coordinate r, we change the variable f —> 
1/x, obtaining V 2 (f> = l/f 4 (d 2 (j/dx 2 ). Substituting this 
in Eq. C9> we have 


d 2 4> dWd 2 W 
dX' 2 ~ dj> d4> 2 ’ 

a 1-dim equation on the x variable. Note that the solu¬ 
tions of the first-order differential equation 


4 = dW 

dx d(j) 

are also solutions of the equation of motion. 

The class of models 

= = ( 5 ) 

is described by the parameter p which drives the way 
the field self-interacts. For p = 1 we have the usual <(> 4 
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theory. For p = 3,5,... the models describe potentials 
supporting minima at <p = 0 and ± 1 , and the solutions 
of Eq.QJ connect the minima <p = ±1 crossing (p = 0, 
on the form of a 2-kink on the x variable. It was shown 
in 0 that the p = 4,6,... cases, where the potential 
has minima at <p = 0 and 1 , lead to stable solutions for 
Eq.®, valid for </> S [0, oo), in the form 


(p(f) = tanh p 


1 

pr 


( 6 ) 


These solutions with even p, p > 4, have spherical sym¬ 
metry, with a central core growing with p. Here we are in¬ 
terested on analyzing the gravitational properties of one 
of these solutions with spherical symmetry. We fix from 
here on p = 4, whose solution on the flat space given 
by Eq. ©, leads to a smooth interpolation between the 
minima <p = 1 and <p = 0 for f going from 0 to oo, re¬ 
spectively. 



FIG. 1: Total energy at a distance R for a defect p = 4 in flat 
space, as a function of R. 

For this case we calculate the energy of the configura¬ 
tion at a distance R from the defect. The result, obtained 
for flat space and depicted in FiglU shows that the energy 
of the defect grows until a plateu is achieved. In this way 
there is no the characteristic divergence obtained for the 
global monopole. This shows that this defect, contrary 
to the global monopole, does not need a distance cutoff 
in order to achieve finite energy. 

As this defect was obtained in flat space, a next step 
is to analyze its gravitational field. With this aim in Sec. 
II we obtain the Einstein and field equations. We show 
that the system of equations need four boundary condi¬ 
tions, two of them related to the behaviour of the scalar 
potential at the center of the core and at the infinity and 
two other ones coming from imposing a Minkowski space- 
time far from the defect. The possible ways to solve com¬ 
pletely this system of equations are discussed. In Sec. Ill 
we consider the Newtonian approximation which shows 
the repulsive character of the defect. We also present an 
analysis of the weak field approximation, leading as a re¬ 
sult that the metric inside the core of the defect presents 


a deficit solid angle. Our conclusions are presented in 
Sec.IV. 


II. EINSTEIN AND FIELD EQUATIONS 


In curved spacetime the action for the defect structure 
is 


S = 


d A 3 


Va 


\ R+ 
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with the usual relation between the spherical coordinates 
r,0,(f> and the Cartesian coordinates ad, i = 1,2,3 and 
where g = — det (g^) and /i, v = 0,1, 2, 3. The potential 
is 

fw = a y7i^ a/2 <’' I,! --A 1/2)2 . (8) 

which is the potential from Eq. © for p — 4, after restor¬ 
ing the dimensional quantity 77 , the vacuum expectation 
value (v.e.v) of the scalar field cp, and the coupling con¬ 
stant A. 

As an Ansatz, we try a general static metric with 
spherical symmetry: 

ds 2 = B(r)dt 2 — A(r)dr 2 — r 2 (dd 2 + sin 2 0d(p 2 ). (9) 

Defining (p = (p/p and f = r/5 , with 8 = A -1 / 2 ?? -1 , the 
equation of motion for the scalar field can be written as 


d 2 cp ( 1 dB 1 dA 2\d(p 
dr 2 \ 2 B dr 2 A dr + r) dr 

- A( 1 _0l/2 ) (30l/2_5^ ) = O (10 ) 

The energy-momentum tensor of the defect is related to 
the variation of the Lagrangian density of the matter field 
with respect to the metric Il5l : 


T = 


2 dC mo 

V9 dg^ v 


( 11 ) 


Using the result dg/dg^ u = — gg (v. 0 , p.364), we 
have for a Lagrangian with an usual (1/2 )d /J ,<p a d IJ, fp a ki¬ 
netic term: 


= d^ a d„r - 9^:1 (12) 

where Cl^ t is the usual Lagrangian of matter field in 
flat space. This gives for our model the following non¬ 
null components for the energy-momentum tensor: 

n = T e= T t = ^ + ^ V ^ ( 13 ) 

T; = -JA^ + ^ VW ' (14) 
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where the prime denotes derivation with respect to r. 
The Einstein equations are 

K - \ R< Jl = g ^ = 8-gt;, (is) 

where 17, and G v are, respectively, the Riemann and 
Einstein tensors. 

The components of the Einstein tensor are easily ob¬ 
tained using the computer algebra systems GRTensorll 
(for Maple): 


G\ 

G r 
r 


A' 1 1 

B' 1 1 

ABr r 2 Ar 2 ’ 

= M = _ 

* 2 ABr 2 A 2 r 

B' 2 


4 AB 2 ' 


(16) 

(17) 

B'A' B" 

AA 2 B + 2 AB ( ' 18 ' 1 

(19) 


With the components of Einstein and energy- 
momentum tensors, and changing to the f variable we 
get explicitly: 


dA 

dr 


dB 

dr 


= A 2 f 


1 

Ar 2 


+ A 


2 A \ dr 
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2 r 4 
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dA B _ / dcj) 

= —pr-Z + ABr ( T- 
dr A \ ar 
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( 21 ) 


and 
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2 ABr 
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2 AB AAB 2 


+ ^ 3/2 (W 1/2 ) 2 )(22) 


with A = r] 2 /M 2 . We see that the Einstein equations 
lEas.CT-d^l together with the equation of motion 
iEq.mil for the <p field form a highly nonlinear set of 
equations. 

The solution <fi(f) for flat space (Eq.©, with p = 4) 
impose the boundary conditions (b.c.) </>( 0 ) = 1 and 

<j>(oo) = 0. Also, as the matter content (related with 
the <p field) decreases with distance, we expect that far 
from the defect the spacetime tend to be four-dimensional 
Minkowski, M 4 . This leads to the additional conditions 
A(oo) = l,R(oo) = 1 . 

From the b.c. we see that the set of equations for A, B , 
c t> is over-determined, and one of the Einstein equations 
is used as a consistency check. In a numerical search 
for a complete set of solutions the b.c. favor the use of 
Ens. ffll - ^rril together with Eq. m- The hybrid nature 
of the b.c. signals that the use of relaxation (or shooting) 
method for Eq. itTUI) together with a Runge-Kutta method 


for Eqs. CT-isn can be used in a search for a solution. 
This study is in progress. In the present work, however, 
as we are interested in a weak field limit, we will use 
Eqs.J^ni and PI together with the solution for the </> 
field obtained in flat space (Eq. 0 for p = 4). 


III. LINEARIZED GRAVITY 


Consider first the Newtonian approximation. The 
Newtonian gravitational potential $(r), is obtained as 
a solution of the Poisson equation 

V 2 < 1? = 8 ttG(T£ - i^T) = 47 tG(Tq - if). (23) 


Substituting the stress-energy tensor on the Poisson 
equation, after setting A = 1 for the Newtonian approx¬ 
imation, we get 


V 2 $ = - 


M 2 r A 


(24) 


or, in terms of the dimensionless variables <j) and f, 


d 2 $ 

'dF 2 


2 0 3 / 2 

^ = -aL-i-. 

r 2 r 4 


ri/2x2 


r- 


(25) 


In this approximation (which for our problem means 
f >> 1 ), we can consider 4>{f) as the solution on flat 
space, 4>(f) = tanh 4 (l/(4r)). A graphic of the r.h.s. of 
Eo. ||251) as a function of r indicates that the defect acts 
as a thick spherical shell with negative mass. In this way, 
we expect a potential with a constant value on the core 
of the defect, with a maximum repulsive force near to 
the shell and decreasing with the distance to the shell. 
As the potential is defined unless a constant, we numeri¬ 
cally solved Eq. PQ) using the second order Runge-Kutta 
method, with the initial conditions $( 0 ) = le $'( 0 ) = 0 . 
The corresponding gravitational force under a particle 
with unitary mass is depicted in Fig HI This figure shows 
a localized aspect of the gravitational force, with the 
magnitude increasing with the parameter A. 

The next step is to investigate the influence of this de¬ 
fect over the spacetime geometry, on the linear approx¬ 
imation. With this purpose we consider A and B from 
Eq.® in the form 


A(r) = 1 + ea(r ), B(r) = 1 + e/3(r). (26) 


The 99 component of the Einstein equations gives, for 
this approximation, 

i^ + ?/5=-^(l-AV. (27) 

After comparison with the Poisson equation for the New¬ 
tonian gravitational potential, we are led to (3 = 24> + c, 
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FIG. 2: Repulsive Newtonian force F(r/5) acting on a particle 
with unitary mass for A = 0.5, 1.0 and 2.0. Here and in 
the following figures the thickness of the lines increases with 
increasing A. 


where c is a constant determined with the boundary con¬ 
dition B( oo) = 1. The result is depicted in Fig|3Ja). 
Also, Eq. <EJ in this approximation leads to 


^ = _^ + A f A 2 

dr dr 'dr 


(28) 


We numerically integrate this equation with <f> given by 
the solution on flat space, and with the result for (3{r) ob¬ 
tained previously. In this way we obtain FigHb). Note 
that Figs|3a)-(b) show that near to the origin, the solu¬ 
tions for A and B tend to different values, indicating the 
appearance of a deficit solid angle. 



FIG. 3: Metric coefficients (a) _B(r/<5) (up) and 

(b)A(r/<5)(down) , on the linear approximation, for A = 0.5, 
1.0 and 2.0. 


IV. CONCLUSION 

In this work we analyzed the gravitational field pro¬ 
duced by a global defect with an explicit 1 /r 4 depen¬ 
dence on the scalar potential. We found that with such 
a potential the defect produced has finite energy. This 
is an alternative to the global monopole 0 0 ) for con¬ 
structing a (3, l)-dim defect, where no cutoff distance is 
needed. The analysis of the Newtonian potential shows 
no effect besides a repulsive gravitational field appear¬ 
ing near the core of the defect. The repulsive Newto¬ 
nian force is proportional to the v.e.v. 77 of the scalar 
field. The analysis of the gravitational field on the lin¬ 
ear regime shows a space with a deficit solid angle. The 
amount of this deficit is controlled by the v.e.v. of the 
scalar field. In fact, Figs|3{a)-(b) show that inside the 
defect, as the parameter A = 77 /Mp increases, the met¬ 
ric coefficients A and B are more departed from unity. 
Our weak field analysis shows the 1/r 4 dependence of 
the scalar potential fixes qualitatively the gravitational 
effects of the defect. The magnitude of these effects is 


governed mainly by the v.e.v. of the scalar potential. 
We expect that a further numerical analysis of the com¬ 
plete set of Einstein and field equations can reveal other 
aspects not presented in this work. 
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